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Hydrodynamic Simulations of Merging Galaxy Clusters: 
Non-Equilibrium Ionization State and Two-Temperature Structure 
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Abstract 

We investigate a non-equilibrium ionization state and an electron-ion two-temperature structure of the 
intracluster medium (ICM) in merging galaxy clusters using a series of N-body and hydrodynamic simula- 
tions. Mergers with various sets of mass ratios and impact parameters are systematically investigated, and 
it is found that, in most cases, ICM significantly departs from the ionization equilibrium state at the shock 
layers with a Mach number of 5-2.0 in the outskirts of the clusters, and the shock layers with a Mach 
number of ^2-4 in front of the ICM cores. Accordingly, the intensity ratio between FeXXV and FeXXVl 
Ka line emissions is significantly altered from that in the ionization equilibrium state. If the effect of the 
two-temperature structure of ICM is incorporated, the electron temperature is ~10-20 % and ~30-50 % 
lower than the mean temperature of ICM at the shock layers in the outskirts and in front of the ICM cores, 
respectively, and the deviation from the ionization equilibrium state becomes larger. We also address the 
dependence of the intensity ratio on the viewing angle with respect to the merging plane. 

Key words: X-rays: galaxies: intergalactic medium — X-rays: galaxies: clusters 



1. Introduction 

Galaxy clusters host an X-ray emitting hot gas with 
temperature of ~ 10 7-8 K, or the intracluster medium 
(ICM), as well as dark matter and galaxies. It is theo- 
retically believed that ICM has been heated up through 
various magnitude of shock waves caused by successive 
mergers of galaxies, galaxy groups and galaxy clusters 
during the hierarchical structure formation (e.g., Ryu et 
al. 2003). These shocks are regarded as sources of tur- 
bulent motion of ICM and sites of particle acceleration. 
Furthermore, they amplify the intracluster magnetic field 
and can be sources of non-thermal emission. In this sense, 
merging galaxy clusters are of great importance as fruit- 
ful sites of interesting phenomena of astrophysical plasma 
(see Sarazin 2002 for a review). 

From an observational point of view, merging galaxy 
clusters provide lots of information about outskirts of 
galaxy clusters, because compression and shock heating 
of ICM during merging events increase the X-ray lumi- 
nosity of ICM at the outer parts of galaxy clusters. For 
example, based on the Suzaku XIS observation, Fujita et 
al. (2008) reported the detection of Fe K emission lines in 
the linked region of Abell 399 and Abell 401, more than 
~ 1 Mpc apart from the centers of both galaxy clusters, 
and estimated the metallicity relatively higher than con- 
ventionally expected. Their detection is in virtue of the 
fact that ICM in the linked region is significantly com- 
pressed, as well as the high sensitivity of the XIS detec- 
tor. Therefore, merging galaxy clusters can be good sites 
for elucidating physical properties of ICM, such as the 



temperature and metallicity, especially in the peripheral 
regions of galaxy clusters. 

Measurements of physical properties of ICM, such as 
the temperature and metallicity, in X-ray observations, 
however, usually assume that ICM is in the ionization 
equilibrium state, and that electrons and ions share the 
same thermal temperature. It is widely believed that these 
assumptions are reasonable inside galaxy clusters, because 
the timescale for collisional ionization equilibrium, 
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10 3 cm -3 

(for e.g., FeXXV and Fexxvi in a temperature of several 
keV; see also Masai 1984) and the timescale for thermal 
equilibration between electrons and ions through Coulomb 
scattering, 



5 x 10 7 K 



3/2 



(2) 



(Spitzer 1962) are usually shorter than the typical dynam- 
ical timescale around the central regions of galaxy clusters, 
where n e and T e are the number density and tempera- 
ture of electrons, respectively. But these timescales are 
not short enough compared with the merger timescale, 
and even longer in the outskirts of galaxy clusters where 
n c ~ 10~ 4 cm~ 3 . Consequently, the assumptions of the 
ionization equilibrium and the thermal equipartition be- 
tween electrons and ions are no longer justified. Therefore, 
in measuring the temperature and/or metallicity of ICM 
around shock layers and outskirts of merging galaxy clus- 
ters, the deviations from the ionization equilibrium and 
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the thermal equipartition have to be taken into consider- 
ation (see also Prokhorov 2009). 

Previously, many works were devoted to numerical 
simulations of merging galaxy clusters under idealized 
but well-controlled initial conditions (Roettiger, Lokcn, 
Burns 1997; Packer, Sarazin 2001; Ritchie, Thomas 2002; 
Takizawa 2005; Takizawa 2006). But most of them ig- 
nored the effects of the non-equilibrium ionization state 
and the deviation from thermal equipartition between 
electrons and ions, or the two-temperature structure of 
ICM. The former effect in the Warm-Hot Intergalactic 
Medium (WHIM) is addressed in cosmological simula- 
tions for the first time by Yoshikawa, Sasaki (2006) and 
Cen, Fang (2006). Takizawa (1999) and Takizawa (2000) 
took into account the latter effect in numerical simula- 
tions of merging galaxy clusters for the first time, while 
Yoshida et al. (2005) and Rudd, Nagai (2009) investigated 
it in WHIM and in outskirts of galaxy clusters, respec- 
tively, in cosmological simulations. Recently, Akahori, 
Yoshikawa (2008) (AY08) carried out an N-body/SPH sim- 
ulation of an Abcll 399 and Abcll 401 system incorporat- 
ing these two effects simultaneously and self-consistently, 
and showed that the non-equilibrium ionization state and 
the two-temperature structure of ICM take place along 
shock layers at the merging interface. 

While AY08 addressed only a particular situation of 
Abcll 399 and Abell 401, in this paper, we systematically 
investigate the non-equilibrium ionization state and the 
two-temperature structure of ICM for various situations 
of merging galaxy clusters, e.g., mass ratios, impact pa- 
rameters, and merging stages, so as to be confronted with 
future X-ray observations of actual merging galaxy clus- 
ters. In addition to that, we also address the viewing 
angle dependence of observational signatures of the de- 
viations from the ionization equilibrium and the thermal 
equipartition. 

The rest of this paper is organized as follows. In section 
2, we describe the numerical methods and the initial con- 
ditions of the simulations. The simulation results are pre- 
sented in section 3. In section 4, we summarize our conclu- 
sions and discuss future prospects for the detection of the 
non-equilibrium ionization state and the two-temperature 
structure of ICM. In the case that the cosmological scaling 
is required, we assume the density parameter £7m = 0.24, 
the baryon density parameter fib = 0.04, and the Hubble 
constant Hq = 70 km s _1 Mpc -1 , unless otherwise speci- 
fied. 

2. Simulation 

2.1. Numerical Method 

We carry out a set of N-body and SPH simulations 
of merging galaxy clusters for various initial conditions. 
In all simulations, radiative cooling is ignored, since, in 
this work, we focus on relatively outer regions rather 
than the very central regions of galaxy clusters, and 
the actual cooling timescale in the central regions of 
the simulated ICM is ~ 5 Gyr, sufficiently longer than 
the merging timescale. Simulations are performed us- 



ing the code developed in the previous papers (Akahori, 
Masai 2005, Akahori, Masai 2006, and AY08), which 
adopts the entropy conservative formulation of SPH by 
Springel, Hcrnquist (2002), and the standard Monaghan- 
Gingold artificial viscosity (Monaghan, Gingold 1983) 
with Balsara limiter (Balsara 1995). In addition to the 
standard SPH prescription, the time evolution of the two- 
temperature structure and the non-equilibrium ionization 
state of ICM are incorporated. 

For the time evolution of the two-temperature struc- 
ture, only the Coulomb scattering is considered as a phys- 
ical process responsible for the thermal relaxation be- 
tween electrons and ions. In this picture, timescales on 
which each of electrons and ions achieves the thermal re- 
laxation are much shorter than that between electrons 
and ions. Therefore, electrons and ions can have differ- 
ent temperatures, T e and Ij, respectively, just after expe- 
riencing shock heating, because nearly all of the kinetic 
energy of ICM is carried by ions at a shock front, and is 
preferentially converted to the thermal energy of ions in 
the post-shock regions. Afterward, thermal relaxation be- 
tween electrons and ions proceeds through Coulomb scat- 
tering on a timescale described by 

3/2 
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l8 vr (r c /io 8 K) 3 / 2 /40_ 

Vr (n;/10- 3 cm- 3 ) ' VlnA 

where m e and vn\ are the electron and ion masses, respec- 
tively, rii the number density of ions, the Boltzmann 
constant, and InA the Coulomb logarithm (Spitzer 1962). 

Assuming that electrons and ions share the same fluid 
velocity, and that shock heating are exclusively con- 
tributed to the heating of ions, the energy equations for 
electrons and ions are given as 



p—7- = -P c V • v + Qex, 
at 

dui 

p—rr = -P\V ■ v + Q sh - Q ox , 
at 



(4) 
(5) 



respectively, where p is the gas mass density, u the specific 
thermal energy, P the thermal pressure (quantities with 
subscript "e" and "i" denote those of electrons and ions, 
respectively) , Q s h the shock heating rate per unit volume, 
and <5 CX the energy exchange rate per unit volume through 
Coulomb scattering: 



Qex — 



n e k Tj-T e 

7-1 *ei 



(6) 



where 7 = 5/3 is the adiabatic index. Using equations 
(4) and (5), and assuming the constant mean molecular 
weight, /i = 0.6, for almost fully ionized ICM, we obtain an 
ordinary energy equation for mixed fluid with a mean tem- 
perature, T = (n c T c + n\T\)/(n c + m), and another equa- 
tion for the time evolution of electron temperature, which 
can be written as 



dt 



Ti - T c T c 



Qsh, 



(7) 
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where T e = T e /T and T\ = T\jT are the dimensionless tem- 
peratures of electrons and ions normalized by the mean 
temperature, respectively, and u — u c + u\. This equa- 
tion is numerically solved in the same manner described 
in Takizawa (1999) as follows. We first regard the sec- 
ond term on the right hand side as being negligibly small, 
and integrate the first term implicitly using the method 
described in Fox, Loeb (1997). Then we integrate the sec- 
ond term explicitly, and finally add this contribution to 
the time evolution of T c . In each integration, the time 
step regulated by the courant condition of SPH ensures 
sufficiently smaller Q s h/u than unity and thus enable us 
to solve the evolution of T e with reasonable accuracy. 

Actually, a possibility is argued that plasma waves are 
able to attain the thermal relaxation between electrons 
and ions more quickly than the Coulomb collision. In 
reality, it is observationally suggested that the electron 
and ion temperatures are almost the same around the 
shock front with a Mach number of ~ 10 — 20 in supernova 
remnants, indicating that physical processes other than 
Coulomb collision contribute to the quick thermal relax- 
ation between electrons and ions (Ghavamian ct al. 2007). 
There exists, however, a room for a very slow shock with a 
Mach number of M < 5 to have significant difference in the 
temperatures between electrons and ions. Thus, in order 
to bracket the plausible range of theoretical uncertainties, 
we present the results both for the two-temperature runs 
described above and the single-temperature runs, in which 
thermal energies of electrons and ions are assumed to be 
instantaneously equilibrated, and electrons and ions share 
a single temperature. 

The time evolution of ionization fractions of ions is com- 
puted by solving the rate equation, 

df- z+1 

k=l 

-ajfj + cxj+ifj+t, (8) 

for each SPH particle, where j is the index of a particular 
ionization stage considered, Z the atomic number, fj the 
ionization fraction of an ion j , Sij the ionization rate of an 
ion j with the ejection of i electrons, and ay is the recom- 
bination rate of an ion j. The ionization processes include 
collisional, Auger, charge-transfer, and photo-ionizations, 
and recombination processes are composed of radiative 
and dielectronic recombinations. The ionization and re- 
combination rates are calculated by utilizing the SPEX 
vcr 1.10 software package 1 . 

Equation (8) is integrated implicitly in the same manner 
described in Yoshikawa, Sasaki (2006), in which ionization 
fractions of all ionic species are sequentially updated in the 
order of increasing ionization states rather than being up- 
dated simultaneously, and the source terms contributed 
by the ionization from and recombination to the lower 
states are evaluated at the advanced time step. In the two- 
tcmpcraturc runs, the reaction rates are computed using 
the electron temperature, T e , calculated from equation (7) 

1 http:/ /www. sron.nl/divisions/hea/spex/ 



Table 1. Simulation runs and their merging parameters 



run ID 


M l [M ] 


M 2 [M ] 


V{ [km s x ] 


b [Mpc 


1:1 head on 


8 x 10 14 


8 x 10 14 


667 


0.0 


1:1 offset 


8 x 10 14 


8 x 10 14 


675 


0.955 


4:1 head on 


8 x 10 14 


2 x 10 14 


645 


0.0 


4:1 offset 


8 x 10 14 


2 x 10 14 


654 


0.955 



so as to incorporate the effect of the two-temperature 
structure. As for the single-temperature runs, the reaction 
rates are computed using the mean temperature, T . 

We solve the time evolution of the ionization fractions of 
H, He, C, N, O, Ne, Mg, Si, S, and Fe; total 112 ionization 
states for each SPH particle. All the ionization states are 
considered in calculations of the X-ray spectrum and the 
surface brightness in this paper. But in the following, 
we mainly focus on the ionization fraction and the line 
intensities of Fe, because the line emissions of Fe are the 
most prominent in the typical temperature range of ICM 
~ 10 7 - 8 K. 

2.2. Initial Condition 

The initial conditions are setup as follows. We adopt 
the NFW density profile (Navarro, Frenk, White 1997) 
for a dark matter halo, and the /3-modcl density profile 
(Cavaliere, Fusco-Femiano 1976) with (3 = 2/3 for a ICM 
component. The scale radius of the NFW profile is set 
to r s = r2oo/5.16, and the core radius of the /3-model is 
t c = r s/2, where r 2 oo is the virial radius defined as the 
radius within which the mean cluster mass density is 200 
times the present cosmic critical density. The radial pro- 
files of velocity dispersion of dark matter and temperature 
of ICM are computed using Jeans equation and the as- 
sumption of hydrostatic equilibrium, respectively. The gas 
and dark matter particles are distributed out to 1.4 r 2 oo, 
which enables us to follow the dynamical evolution of ICM 
at outskirts of galaxy clusters. 

In this paper, mass ratios of two galaxy clusters are 
set to 1:1 and 4:1. For runs with 1:1 mass ratio, the 
virial mass and radius are set to M 2 oo = 8x 10 14 M Q and 
^200 = 1-91 Mpc, respectively. For runs with 4:1 mass ra- 
tio, a main (or more massive) cluster is the same as the 
one in the runs with 1:1 mass ratio, and the virial mass 
and radius of a sub (or less massive) cluster are set to 
M 20 o = 2 x 10 14 M & and r 2 oo = 1-21 Mpc, respectively. 
Density profiles of dark matter and ICM and temperature 
profiles for the two masses of galaxy clusters are depicted 
in figure 1. For each mass ratio, we carry out the runs 
with different impact parameters, b; "head-on" runs with 
6 = and "off-set" runs with b = 0.955 Mpc. As for the 
metallicity, a spatially uniform mctallicity of 0.2 times the 
solar abundance is assumed. In addition to that, we also 
assume that electrons and ions share the same tempera- 
ture, and that they are in the ionization equilibrium state 
at the start of the simulations. 
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Fig. 1. Radial density profiles of ICM (baryon) and dark 
matter normalized by the cosmic mean density of each com- 
ponent (lower panel), and ICM temperature (upper panel) at 
initial conditions. Solid lines indicate the profiles of galaxy 
clusters in the runs with 1:1 mass ratio and the main cluster 
in the runs with 4:1 mass ratio (M200 = 8 X 10 14 Mq), while 
dashed lines show the profiles of the sub-cluster in the runs 
with 4:1 mass ratio (M200 = 2 X 1O 14 M0). The vertical lines 
indicate the locations of their virial radii. The central ICM 
density corresponds to n a = 5.52 X 10 — 3 cm -3 . 



The initial separation of two galaxy clusters, d, is set 
to 1.4 times of the sum of their virial radii so that they 
contact each other at their outer edge. The initial relative 
velocity of the two clusters, v u is calculated using the 
prescription described in Sarazin (2002), and is given by 



2G(M 1 + M 2 ) 



1 

d 



1 - 



(9) 



where G is the gravitational constant, and Mi and M2 
are the virial masses of the two clusters, do denotes the 
maximum separation where their relative radial velocity 
is zero, and is calculated with the aid of the Kepler's third 
law as follow 

d = [2G(M 1 + M 2 )} 1 / 3 ( t -^) 1 , (10) 



where to is the comic time when the merging event takes 
place, and is set to to = 13.7 Gyr in this paper. All the 
runs presented in this paper and their merging parameters 
are listed in table 1. 

In all simulations, the mass of each dark matter and 
SPH particle is set to 8.62 x 10 8 M Q and 1.64 x 10 8 M Q , 
respectively. The single cluster in the runs with 1:1 mass 
ratio, and the main cluster in the runs with 4:1 mass ratio 
are composed of one million dark matter particles and the 
same number of SPH particles, while the sub-clusters in 
the runs with 4:1 mass ratio are represented by a quarter 
million particles for each component. The correspond- 
ing spatial resolution (the smoothing length) of SPH is 
~ 50 kpc at the ICM density of 10 -3 cm -3 (central part of 



the clusters), or 
clusters) . 

3. Result 



110 kpc at 10 4 cm 3 (outskirts of the 



3.1. Single- Temperature Runs 

First, we show the results for the single-temperature 
runs, in which thermal relaxation between electrons and 
ions is supposed to be achieved very quickly, and electrons 
and ions share a single temperature. 

Figure 2 shows the time evolution of the ICM den- 
sity (color), the dark matter density (black contour), 
and the Mach number distribution of ICM (white con- 
tour), on a collision plane of the two merging galaxy 
clusters, where the Mach number is estimated with the 
shock-capture scheme described in Pfrommer et al. (2006). 
Panels from top to bottom show the snapshots at a time 
of t = -0.5 Gyr, 0.25 Gyr, 0.375 Gyr and 0.5 Gyr, re- 
spectively, where t = Gyr corresponds to the time of the 
closest approach of the centers of the dark matter halos. 
Therefore, dark matter halos at t > have already passed 
each other. 

In the early stage of the merging (£ = — 0.5 Gyr), the out- 
skirts of the two galaxy clusters interact with each other, 
and then form the linked regions where ICM is compressed 
due to the collision. At the edges of the linked regions, 
there exist shock layers with a Mach number of 1.5-2.0 
(white contour in the top panels of figure 2). These shocks 
are immediately generated at the contact interface of the 
two clusters when the clusters start to collide, then prop- 
agate to outside of each cluster. 

The emergence of these shocks at the outskirts of galaxy 
clusters is also pointed out by AY08, but not clearly 
seen in the previous SPH simulations (Takizawa 1999; 
Takizawa 2000). This is because, in these previous stud- 
ies, density profiles of ICM are truncated in the outside 
of clusters (r > 5r c ), and also because these studies adopt 
isothermal temperature profiles of ICM at the start of the 
simulations. On the contrary, in the present study, the dis- 
tributions of ICM are extended out to r = 14.4 r c , and the 
ICM temperature decreases toward outside as suggested 
by recent X-ray observations (e.g., Vikhlinin ct al. 2005). 
Therefore, the sound speed at the outskirts of galaxy clus- 
ters, which is typically ~ 500 km s _1 for a gas tempera- 
ture of ~ 2 keV, is lower than the relative velocity (~700- 
800 km s~ x ) of the two clusters at this stage, resulting in 
the formation of the shock waves. 

At the middle stage of the merging (t= 0.25-0.375 Gyr), 
another shock waves with a Mach number of 1.5-2.0 in 
front of the ICM density peaks emerge, because the dense 
cores of ICM are accelerated by the preceding dark matter 
halos, and they have Mach numbers of ~3-4 and 2-3 at 
the late stage of the merging (t = 0.5 Gyr) in the case 
of 1:1 and 4:1 mass ratios, respectively. Since the shock 
waves formed in front of the ICM cores are faster than the 
ones formed in the outskirts at the early stage, the former 
overtake and merge with the latter, as can be seen in the 
bottom panels of figure 2. 

As shown in the top panels of figure 3, most of ICM ex- 
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Fig. 2. Density maps of ICM on a collision plane of the two galaxy clusters are shown for 1:1 head on, 1:1 offset, 4:1 head on, and 
4:1 offset (from the left to the right) at a time of t = —0.5 Gyr, 0.25 Gyr, 0.375 Gyr, and 0.5 Gyr (from the top to the bottom), 
respectively, where t = Gyr corresponds to the time of the closest approach of the centers of the dark matter halos. White contours 
indicate the Mach number distribution of ICM for a Mach number of 1.4 and above, and are separated by a difference of 0.6. The 
innermost contours for 1:1 and 4:1 mergers at a time of t = 0.5 Gyr are 3.2 and 2.6, respectively. The distributions of dark matter 
are also overlaid with black contours. 
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Fig. 3. The mean temperature of ICM, f, the ratio of the ionization temperature relative to the mean temperature, Tj on /T , and 
the ratio of the ionization fraction of Fe XXV relative to that in the ionization equilibrium state (from the top to the bottom) on a 
collision plane of the two merging galaxy clusters in the runs, 1:1 head on, 1:1 offset, 4:1 head on, and 4:1 offset (from the left to 
the right), at a time of t = 0.5 Gyr. 
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Fig. 4. Ionization fractions of Fe in the ionization equilib- 
rium state as a function of temperature. 

cept the regions behind the shocks is in the temperature 
range of 5 ^ T < 15 keV, and in this temperature range, 
Fexxv and Fexxvi arc dominant among all the ioniza- 
tion states of iron, if ICM is in the ionization equilibrium 
state (see figure 4). Therefore, in the rest of this work, 
we mainly focus on the ionization fractions of Fe xxv and 
Fexxvi to address the degree of the deviation from the 
ionization equilibrium in ICM. In addition, we introduce 
the "ionization temperature", T lon , the temperature at 



which given ionization fractions of Fe xxv and Fe XXVI are 
realized in the ionization equilibrium state to characterize 
the ionization state of ICM. Note that T; on monotonically 
increases with the ratio of abundances of Fe xxvi relative 
to Fexxv, and differs from T if these ions are not in the 
ionization equilibrium state. Thus, the difference between 
2} on and T can be used as a measure of the deviation from 
the ionization equilibrium state. Considering the situation 
where ICM is shock-heated from a few keV to higher tem- 
perature, Tj on is lower than T just after the shock heating 
since the ionization from Fexxv to Fexxvi is not quick 
enough to catch up with the ionization equilibrium. After 
that, Tj on asymptotically approaches to T as the ioniza- 
tion state gets close to the ionization equilibrium. 

The middle and bottom panels of figure 3 show the 
maps of the ratio of the ionization temperature relative 
to the mean temperature of ICM, Tj on /T, and the ra- 
tio of the fractions of Fe XXV relative to that in the ion- 
ization equilibrium state, /(Fexxv)// oq (Fexxv), at the 
late stage (t = 0.5 Gyr), respectively. Note that the frac- 
tion in the ionization equilibrium, / eq (Fexxv), is calcu- 
lated using the mean temperature, T, It is clearly seen 
that we have a significant difference between T- lon and T 
around the shock layers depicted in figure 2, indicating 



No. ] 

that there exists a significant deviation from the ioniza- 
tion equilibrium state. The ratio of the ionization frac- 
tions, / (Fe xxv) // Gq (Fexxv), is above unity and ~ 1.5 at 
the shocks in the outskirts, and ^5 — 7 and ~ 4 at the 
shocks in front of the ICM cores in the runs with 1:1 and 
4:1 mass ratios, respectively. 

The over-population of Fe xxv fraction at the shocks in 
the outskirts and in front of the ICM cores described above 
can be understood as follows. In these regions, the shock 
heating increases the ICM temperature from 2 — 3 keV to 
the higher temperature, so that the Fe XXV fraction in the 
post-shock regions declines with time toward the equilib- 
rium fraction because the Fe xxv fraction is maximum at 
~ 3 keV in the ionization equilibrium state (see figure 4) . 
However, the ionization of Fexxv toward higher-ionized 
states is not quick enough to catch up with the ionization 
equilibrium state. Therefore, the Fexxv fraction is higher 
than that in the ionization equilibrium state. 

3.2. Two-Temperature Runs 

In this subsection, we show the results for the two- 
temperature runs, in which the assumption that electrons 
and ions share a single temperature is relaxed, in order 
to examine the effect of the two-temperature structure on 
the ionization state of ICM. The thermal relaxation be- 
tween electrons and ions is solved using equation (7). The 
simulations are performed for the same set of the initial 
conditions as the single-temperature runs. 

The upper panels of figure 5 depict the maps of the ratio 
of the electron temperature relative to the mean temper- 
ature of ICM, T c /T, on the collision plane at a time of 
t = 0.5 Gyr. It can be seen that the electron temperature 
is typically <~ 10-20 % lower than the mean temperature 
at the shocks in the outskirts, and ~ 50 % and ~ 30 % 
lower at the shocks in front of the ICM cores in the runs 
with a mass ratio of 1:1 and 4:1, respectively. In the runs 
with 1:1 mass ratio, we have a higher Mach number and 
a larger jump of the ion temperature at the shock layers 
than the runs with 4:1 mass ratio, and accordingly the 
larger temperature difference between electrons and ions. 

Since the ionization and recombination rates in equa- 
tion (8) depend on the electron temperature, the ioniza- 
tion states at the shocks in the two-temperature runs 
are altered from the results in the single-temperature 
runs in which the ionization and recombination rates 
are computed using the mean temperature of ICM, T. 
The lower panels of figure 5 show the ratio of Fexxv 
fractions / 2T (Fexxv) in the two-temperature runs rel- 
ative to / eq (Fexxv) computed using T. Comparing 
the bottom panels of figure 3 and 5, we can see that 
the ratios, / 2T (Fe XXV)// eq (Fe xxv), are larger than 
/ (Fe xxv) // oq (Fexxv) in the single-temperature runs at 
the shocks in front of the ICM cores. This is because 
a delay of the heating of electrons at the shocks leaves 
Fe xxv fraction higher than that in the single-temperature 
runs. In other words, the electron temperature in the two- 
temperature runs (10-13 keV) is lower than that in the 
singlc-tcmpcrature runs (15-25 keV), and thus the ion- 
ization rate of Fexxv to the higher ionization levels is 
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significantly smaller. 

3.3. Line Intensity Ratio 

Here, let us consider observational imprints of the de- 
viation from the ionization equilibrium and the difference 
in temperature between electrons and ions in the X-ray 
spectroscopic observations. X-ray spectra of the simulated 
galaxy clusters are computed with the simulated ioniza- 
tion state and the hydrodynamical properties of ICM by 
utilizing the SPEX vcr 1.10 software. In order to char- 
acterize the signatures of the non-equilibrium ionization 
state, we focus on the intensity of Fe Ka line emissions 
which is usually adopted as a diagnostic feature of the 
ICM mctallicity in X-ray observations, and introduce a 
ratio of the X-ray intensity between the two energy bands 
defined as 

7(6.6- 6.7 keV) 

/(6.9-7.0keV)' [ ' 

where 7(6.6 — 6.7 keV) and 7(6.9 — 7.0 keV) are mainly con- 
tributed by the intensities of Ka emissions from Fe xxiv- 
XXV and FeXXVi, respectively. Since the contribution of 
Fexxiv Ka emission to 7(6.6 — 6.7 keV) is minor when 
/(Fexxiv) < 0.2, 7? roughly reflects the ratio of Ka line 
emission intensity between Fexxv and Fexxvi. Then, 
we introduce a ratio, R/ R eq , as a plausible observational 
tracer of the deviation from the ionization equilibrium, 
where 7? oq is the intensity ratio under the assumption of 
the ionization equilibrium and thermal equipartition be- 
tween electrons and ions. Note that R/R cq is independent 
of the local ICM metallicity. 

Figure 6 shows the maps of R/R cq of the simulated 
clusters at a time of t = 0.5 Gyr in the single- and two- 
temperature runs viewed along the line of sight perpen- 
dicular to the collision plane. It is clearly seen that 
R/R cq significantly exceeds unity at the shocks in the 
outskirts and in front of the ICM cores. This result is 
consistent with the fact that Fexxv is over-populated 
compared with the ionization equilibrium state at these 
shocks layers. The deviation of R/R cq from unity in the 
two-temperature runs is larger than that in the single- 
temperature runs in accordance with the fact that the 
departure from the ionization equilibrium state in the 
two-temperature runs is larger than that in the single- 
temperature runs. Furthermore, since the difference in 
temperature between electrons and ions is larger at the 
shocks in front of the ICM cores than at the shocks in 
the outskirts, we have significant difference in R/R cq be- 
tween the single- and two-temperature runs primarily at 
the shocks in front of the ICM cores. 

Maps of R/R cq in the two-temperature runs viewed 
from different viewing angles (30 and 60 degrees) are 
shown in figure 7. Compared with the maps of R/R cq 
for a viewing angle of degree in figure 6, the deviation of 
R/ R cq from unity is significantly smaller for a viewing an- 
gle of 30 degrees, and almost negligible for a viewing angle 
of 60 degrees. This is because the shock layers where the 
ionization state deviates from the equilibrium state have 
geometrically thin structure and are almost perpendicu- 
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Fig. 5. The ratio of the electron temperature relative to the mean temperature of ICM, T c /T, (upper) and the ratio of the ionization 
fractions of Fe XXV relative to that in the ionization equilibrium state (lower) on a collision plane of the two merging galaxy clusters 
in the runs 1:1 head on, 1:1 offset, 4:1 head on, and 4:1 offset (from the left to the right), at the late stage of the mergers (t = 0.5 Gyr) 
in the two-temperature runs. 



1:1 head on 1:1 offset 4:1 head on 4:1 offset 




x [Mpc] x [Mpc] x [Mpc] x [Mpc] 

Fig. 6. The ratio, R/R cq for the runs 1:1 head on, 1:1 offset, 4:1 head on, and 4:1 offset (from the left to the right), at t = 0.5 Gyr 
in the single- (upper) and two-temperature (lower) runs. 
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Fig. 7. Same as figure 6 but with viewing angles of 30 degrees (upper) and 60 degrees (lower) in the two-temperature runs. 



lar to the collision plane. Therefore, when viewed with a 
finite viewing angles, the signatures of the deviation from 
the ionization equilibrium are diluted by the surrounding 
ICM which is almost in the ionization equilibrium state. 

3.4- Projected Temperature and X-ray Surface 
Brightness 

Finally, let us consider the projected temperature and 
X-ray surface brightness in the merging galaxy clusters. 
Observed temperature of ICM estimated by the X-ray 
spectrum fitting is generally different from its physical 
temperature due to the projection effect along the line 
of sight and the characteristic response of X-ray detec- 
tors. Here, we adopt the spcctroscopic-like temperature 
(Mazzotta et al. 2004) as an estimate of the observed pro- 
jected temperature of ICM. 

Figure 8 shows the maps of the spectroscopic-like tem- 
perature and the X-ray surface brightness in the 2-10 keV 
energy band viewed from different viewing angles at a 
time of t = 0.5 Gyr. For a viewing angle of degree, 
jumps of the projected temperature across the shocks are 
~ 7 — 13 keV in front of the ICM cores and ~ 3 — 7 keV in 
the outskirts. When viewed with a viewing angle of 30 or 
60 degrees, the jumps of the projected temperature across 
the shock fronts are blurred and smaller than those in the 
case with a viewing angle of degree, because the shock 
fronts are not aligned with the line of sight. Therefore, the 
Mach numbers estimated from the jumps of the projected 
temperature could be biased. Actually, the Mach num- 
bers inferred from the apparent jumps of the projected 
temperature for a viewing angles of 0, 30, and 60 degrees 
in the case of 1:1 head on run are 3.8, 2.2, and 1.5 while 
the actual Mach number is ~ 4. The jumps of projected 
temperature in the single-temperature runs are almost the 
same as those in the two-temperature runs, because the 



signatures of the two-temperature structure are diluted by 
the surrounding ICM which is in the thermal equiparti- 
tion, when viewed with a finite viewing angle. Therefore, 
we can conclude that these biases in the Mach numbers 
are mainly due to the projection effect. 

As for the X-ray surface brightness, sharp discontinu- 
ities at shock layers are clearly visible only for a viewing 
angle of degree. This is again because the shock fronts 
are geometrically thin and the projection effects dilute 
the discontinuities in the case of non-zero viewing angles. 
On the other hand, contact discontinuity surfaces for 4:1 
runs are visible in X-ray surface brightness when a viewing 
angle is degree and even 30 degrees, since the contact 
discontinuity surfaces for 4:1 runs have "bullet" -like struc- 
tures and their appearance in the X-ray surface brightness 
relatively insensitive to the viewing angles. 

4. Summary and Discussion 

We carry out N-body + SPH simulations of merging 
galaxy clusters with various sets of initial conditions by re- 
laxing the assumptions of both the ionization equilibrium 
and thermal equilibration between electrons and ions in 
order to investigate the ionization state of ICM in the 
merging galaxy clusters in detail. As for the thermal 
equipartition between electrons and ions, we perform the 
single- and two-temperature runs to bracket the plausible 
theoretical uncertainties regarding the thermal relaxation 
processes between electrons and ions. 

In our simulations, wc find that the shocks induced by 
the merging events can be classified into two types regard- 
less of the initial conditions. The first ones are formed at 
the early stage of the merging processes in the outskirts 
of the clusters where the sound velocity is roughly equal 
to or lower than the initial relative velocity of the merging 
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Fig. 8. The spectroscopic-likc temperature map viewed with a viewing angle of 0, 30, and 60 degrees (from the top to the bottom), 
respectively, for the runs 1:1 head on, 1:1 offset, 4:1 head on, and 4:1 offset (from the left to the right), at t = 0.5 Gyr in the 
two-temperature runs. Contours shows the maps of X-ray surface brightness in the 2—10 keV energy band. 



galaxy clusters, and the second ones are formed in front 
of the ICM cores in the middle and late stages. Since the 
cores of ICM is accelerated by the gravitational poten- 
tial of collisionless dark matter, the Mach number of the 
shocks in front of the ICM cores is higher than that in the 
outskirts. 

It is shown that the ionization state of iron around these 
shocks deviates from the ionization equilibrium state. 
More specifically, the fraction of Fexxv, the most pop- 
ulated ionization state of iron in typical galaxy clusters, 
exceeds its ionization fraction in the ionization equilib- 
rium at the shock layers both in the outskirts and in front 
of the ICM cores. This excess of the Fexxv fraction can 
be understood by the facts that the ionization processes 
at the shock layers are not quick enough to achieve the 
ionization equilibrium instantaneously and that the ion- 
ization state there differs from the equilibrium state for a 
while after experiencing the shock heating. Furthermore, 
in the two-temperature runs, the electron temperature at 
the shock layers are remarkably lower than the mean tem- 
perature of ICM, and thus the deviation from the ioniza- 
tion equilibrium state is more significant than that in the 
single-temperature runs, because the ionization processes 
at the shock layers are retarded due to the lower temper- 
ature of electrons. 



From the observational point of view, the deviation 
from the ionization equilibrium can be identified by the 
difference between intensity ratios of Fe Ka line emissions 
and that inferred from the ICM temperature estimated 
by the X-ray continuum spectra. We find that the ratio, 
R/R cq , can be a good tracer of such a deviation from the 
ionization equilibrium. We also find that the appearance 
of such observational imprints of the non-equilibrium ion- 
ization state is very sensitive to the viewing angle of the 
merging galaxy clusters to the observers. R/R eq is the 
largest if the merging galaxy clusters are observed along 
the line of sight perpendicular to the collision plane of the 
two clusters, because the shock layers are geometrically 
thin and are parallel to the line of sight. On the other 
hand, the significance of such imprints is weaker in the 
case of a viewing angle of 30 degrees, and almost negligi- 
ble in the case of 60 degrees. It is found that jumps of the 
spectroscopic-like (projected) temperature and the X-ray 
surface brightness across the shock fronts also depend on 
the viewing angle, and that the Mach numbers estimated 
from the jumps of the observed temperature could be also 
underestimated in the case of a finite viewing angle. 

While we concentrate only on the ionization state of Fe 
so far, detections of the deviation from the ionization equi- 
librium of other heavy elements could provide independent 
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information on physical properties of shock waves. Among 
heavy elements other than Fe, Si is the most interesting 
from the observational point of view because the emis- 
sion lines of Sixni-xiv at a photon energy of ~ 2 keV 
are relatively prominent. Actually, at the shocks in the 
outskirts the ionization fraction of Sixiv is 1.2-1.5 and 
1.2-2.0 times larger than that in the ionization equilib- 
rium state in the single- and two-temperature runs, re- 
spectively. Since the ionization fraction of Si XIV in the 
ionization equilibrium is maximum at a temperature of 
~ 10 7 K, and almost fully ionized at higher temperature, 
the observational signatures of the deviation from the ion- 
ization equilibrium could be seen in the outskirts rather 
than in the central regions of galaxy clusters. 

With the current X-ray observational facilities, the sep- 
arate detection of emission lines of Fe xxv and Fe xxvi, 
or the measurement of the ratio R/R eq is still difficult 
because of a lack of the spectroscopic resolution of the 
CCDs. Therefore, spectroscopic analyses to estimate the 
metallicity of ICM with the assumption that ICM is in 
the ionization equilibrium state can lead to a biased re- 
sult. A good spectroscopic resolution is also crucial to es- 
timate the temperatures of electrons and ions separately, 
because the ion temperature is only estimated from the 
detailed line profiles of emission lines, while the elec- 
tron temperature can be relatively easily estimated from 
the X-ray continuum spectra. Thus, the precision spec- 
troscopy in X-ray observations is of crucial importance 
in studying the non-equilibrium ionization state and the 
two-temperature structure of ICM, and can be achieved 
by X-ray calorimeters on board future satellites such as 
Astro-H and International X-ray Observatory. 
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